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I. INTRODUCTION 

Characterizing the dynamical behavior of complex 
quantum systems is an important yet daunting task. The 
standard approach known as quantum process tomogra- 
phy (QPT) [iHSl can provide full information about the 
quantum process, but it consumes a huge amount of re- 
source. Namely, in order to fully determine a quantum 
operation acting on a system consisting of n qubits, QPT 
needs to use 8(16") observables to estimate all the pa- 
rameters necessary to describe this operation. Even if 
this operation is knovifn to be unitary, it still needs 0(4") 
observables. Although many improvements and variants 
of QPT have been proposed [4-19], in general the re- 
source consumption of this approach still grows quickly 
as the system becomes large. 

On the other hand, in many situations we might not 
need to fully determine the quantum operation, but 
merely wish to know whether it satisfies certain prop- 
erty or is far from having this property (assuming it is 
one of the two cases). For example, given a quantum 
machine acting on an n-qubit system, it is natural to ask 
whether it only nontrivially acts on a few qubits, or it 
has non- negligible effect on every qubit. Similar ques- 
tions have been raised and studied in the classical situa- 
tion. For example, given a boolean function (or a graph) 
as an oracle, we want to know whether the function is 
linear (or the graph is connected) or is far from any of 
such functions (or graphs) with respect to some reason- 
able metric, by making only a few queries to the oracle. 
This problem is usually called property teg^mo [20l - l22j . 
It has been extensively studied in computer science and 
has wide applications such as probabilistically checkable 
proofs (PGP) [13 . Surprisingly, many properties of func- 
tions and graphs are found to be testable with very few 
queries. Sometimes the query complexity is even inde- 
pendent of the input's size. 



Given these facts, one may naturally wonder whether 
the less ambitious goal of property testing of quantum 
states or operations would lead to a dramatic decrease 
in resource consumption. Several previous results indi- 
cated that it is indeed the case |24l - l28l |. For example, 
the separability of multipartite states and operations can 
be tested with query complexity independent of the sys- 
tem's sizepGj. In this paper, we will continue this line of 
research and focus on studying property testing of uni- 
tary operators. (The reader should not confuse our work 
with quantum property testing [29] , which concerns about 
the testing of classical objects with quantum algorithms. 
Here we are interested in the testing of quantum objects 
themselves.) We first introduce a normalized distance 
measure that quantifies the average difference between 
unitary operators. Then we show that, with respect to 
this distance measure, the orthogonal group, quantum 
juntas (i.e. unitary operators that only nontrivially act 
on a few qubits of the system) and Clifford group can 
be all efficiently tested. In fact, their testing algorithms 
have query complexities independent of the system's size 
and also have only one-sided error. Next, we give a gen- 
eral algorithm that tests any finite subset of the unitary 
group, and demonstrate an application of this algorithm 
to the permutation group. This algorithm also has one- 
sided error and polynomial query complexity, but we do 
not know whether it can be efficiently implemented in 
general. 

The remainder of this paper is organized as follows. 
In section 2, we introduce the definitions, notations and 
tools used in this paper. Then, in sections 3, 4 and 5, we 
study the testing of orthogonal group, quantum juntas 
and Clifford group respectively. In section 6, we present 
an algorithm that tests any finite subset of the unitary 
group, and then exhibits its application to the permuta- 
tion group. Finally, section 7 concludes this paper. 
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II. PRELIMINARIES 
A. Definitions and Notations 

In this section, we introduce the definitions and nota- 
tions used in this papers. 

Let n > 1 and N — 2". We use Mat^jv to denote the set 
of hnear operators from to (given a fixed basis 
for C^, they are represented by x iV matrices with 
complex entries), and use U^v = {U € Mjv,jv : UW = 
U^U = 1} to denote the set of A''-dimensional unitary 
operators. We are going to regard M.n,n as a Hilbert 
space equipped with the Hilbert-Schmidt inner product 



{A,B) =tT{A^B) 



(1) 



This inner product induces the Hilbert-Schmidt (or 
Frobenius) norm for A ~ {o-i,j)fj=i 



\A\\ = ^tr(Atyl) 



N N 



1=1 j=i 



This norm further induces the following metric 

d{A,B) ^ \\A-B\\. 



(2) 



(3) 



But this metric might be not good for comparing unitary 
operators, since in general we have d{U, V) ^ d{e^^U, V) 
for 6 £ (0, 27r), although U and e^^U are usually consid- 
ered as the same operation since they are equivalent up to 
a global phase. To overcome this problem, we introduce 
another distance measure as follows. First, wc define an 
equivalence relation between linear operators as follows: 
A ~ B if and only if A = e^^B for some 6 e [0,27r). 
Then for any A £ Mn^n, define [A] = {B : A - B}. 
The distance between A and B is given by 



DiA,B) = 



mm 

ce[A],De[B] 
1 

min 

ee[o,27r) ^2N 




(4) 



More generally, for any S C Mjv,jv, define [S] = UAes[A]. 
And the distance between two sets S and T is given by 

D(5,r) = min D{A,B) 
Aes,BeT 

= min ^=||A-B||. ^ ' 

AG [5], Be [T] ^/2N 

It can be easily checked that 

1. D{A, B) > 0, and the equality holds if and only if 
A^B. 

2. D{A,B) = D{B,A). 

3. D{A, B) + D{B, C) > D{A, C). 

Besides, for unitary operators, D has the following nice 
properties: 



1. D{U,V) < 1. 

2. D{UVi,UV2) ^ D(yi,V2). 

3. D{U(E)Vi,U(E)V2) = D{VuV2). 

4. D{UiVi,U2V2) < D{Uu U2) + D{Vi,V2). This is a 
consequence of D{UiVi,U2V2) < D{UiVi,U2Vi) + 
D{U2Vi,U2V2) and property 2. 

5. D{Ui<S)Vi,U2<»V2) < D{Ui,U2) + D{Vi,V2). This 
is a consequence of D{Ui (g) Vl, C/2 (81 V2) < D{Ui (8) 
Vi,U2® Vi) + D{U2 (81 T4, C/2 (81 V2) and property 3. 

Thus I? is a normalized distance measure that reflects 
the average difference between unitary operators. In ad- 
dition, the following relation between D{U, V) and {U, V) 
would be very useful: 

D'iU,V)^l-^\{U,V)\. (6) 



B. Our Question 

The task of property testing is typically described as 
follows. Suppose some unknown object, such as a graph 
or a boolean function, is given as an oracle which can 
be queried locally many times. Our goal is to determine 
whether this object has certain global property or is far 
from having this property, by making as few queries as 
possible. 

Formally, let i7 be a predetermined set from which the 
object is chosen. f2 should be also equipped with a dis- 
tance measure d. A property is a subset 5 C f2. For 
any A £ H., ii A £ S, then we say that A has property 
S; otherwise, if d{A, S) > e, i.e. d{A, B) > e for any 
B £ S, then we say that A is e-far from property S. An 
algorithm e-tests property S if for any input A £ fl, 

• HA has property S, then the algorithm accepts A 
with probability at least 2/3; 

• if j4 is e-far from property S, then the algorithm 
accepts A with probability at most 1/3. 

Besides, if the algorithm makes at most Q(|i^|, e) queries 
to the oracle, then we say that it has query complexity 
0(g(|f2|, e)). A testing algorithm would be very efficient 
if its query complexity depends only on e but not on \ 

In this paper, we will study the problem of property 
testing of unitary operators. In our case, fl — Vn and 
we use D defined as Eqs.Q and ([5]) as the distance mea- 
sure. However, we need to slightly change the definition 
of having a property and being far from a property as 
follows: let S C Vn be a property. We say U £ Uat 
has property S if D{U,S) = 0, i.e. U £ [5]; otherwise, 
we say U is e-far from property S if D{U,S) > e, i.e. 
D{U, V) > e, VV^ £ S. Our input is a blackbox imple- 
menting some [/ e Ujv which can be accessed as follows: 
first, we prepare some state IV'ab), where A, B are its 
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two subsystems such that dini(A) = N and B is some 
auxihary system; then we apply U on the A subsystem, 
obtain IV'^^) = (C^ ® finahy we perform some 

measurement on IV'ab) ^^'^ S^t information about U. In 
certain cases, we are also allowed to access a blackbox 
implementing W . Our goal is that given any S C Vn, 
find an algorithm that e-tests S with the minimal query 
complexity. 



C. Useful tools 



The following tools will be very useful for our work. 



1. Choi-Jamiolkowskt Isomorphism 

The Choi-Jamiolkowski isomorphism [sO, HJ states 
that there is a duality between quantum channels and 
quantum states. In particular, there exists an isomor- 
phism between unitary operators in Uat and pure states 
[nCNxN_ Specifically, let 



1 N 



(7) 



be the A^-dimensional Bell state. For any A € Mjv,Af, 
define 

\v{A)) = {A®m%), (8) 
where A is applied to the first subsystem. Then we have 



{v{A)\v{B)) ^ j^{A,B). 
In particular, for any U,V € Uat, we have 



1 



{v{U)\v{V)) = j^{U,V), 



(9) 



(10) 



where S = diag(cTi , (72, . . . , (tn) with at > 0, and Vi , V2 G 
Uat. Then we have 



N 



Mil ^ WE.?. 



(13) 



If A is unitary, then its singular values cr^ = 1. The 
following lemma shows that the converse is also true in an 
approximate sense: (In what follows, when we write A < 
B for A,B £ Mtv.at, we means that B — A is semidcfinite 
positive.) 

Lemma 1. If A e Mn,n satisfies A^ A < I and > 
N{1 - e), then D{A, Vn) < 

Proof: Suppose A has the singular value decomposition 
as Ea. (fT2| . Note that the condition A^A<Iis equiva- 
lent to 

a, < 1, Vi. (14) 
And the condition > A^(l — e) is equivalent to 



N 

Eaf>7V(l-e). 

1=1 



(15) 



Define U = ViV2 eUn- Then 



DHU,A) < ^WU-Af 

N 

i=l 
N 



(16) 



< e/2, 

which implies D{U,A) < ^JTJl. 



3. Pauli Decomposition 



and 



Let 



Hu))\\^\\\v{v))\\ 



(11) 



So the angle between \v{U)) and |w(V")) faithfully refiects 
the "angle" between U and V with respect to the Hilbert- 
Schimdt product. And if we perform the projective mea- 
surement {\v{V)){v{V)\,I - \v{V)){v{V)\} on the state 
\v{U)), then the probability of obtaining the outcome 
corresponding to \v{V)){v{V)\ is \{U,V)\'^ /N^ . 



2. Singular Value Decomposition 

Suppose A has the singular value decomposition 

A = ViJ:V2, (12) 



^0 -i 
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^2 = I ^ I = r, a3 = ( J \ \=Z 



(17) 



be the Pauli operators. And let Z4 = {0, 1,2,3}, 
{0, 1, 2, 3}". For any x — (a;i, 2:2, • • • , a;„) G ZJ , define 



(18) 



Then {-^agj^ez^ form an orthonormal basis for M^rjv 
with respect to the Hilbert-Schmidt inner product. So 
we can write any A G Mjv,Ar as 



(19) 
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where 



It can be easily checked that 

\\A\?^N E = 

In particular, any t/ G Uat satisfies 



(20) 



(21) 



Uur- E IMC/)P = 1- (22) 



By the Choi-Jamiolkowski isomorphism, {|f (crj))}:Egzj 
also form an orthonormal basis for C^^^. For any A e 
Mat at, we have 



HA)) = E M.-(^)l«(a;)). 



(23) 



Hence, if we measure the state |f (^)) (assuming it is 
normalized) in the basis {\v{(jg)) : x G ZJ}, we get the 
outcome x with probability |/ij(yl)p. 

Finally, we define ® : Z^f x ^ ZJ and : ZJ x ^ 
Z4 such that for any x,y £ Z" , 



(24) 



Then 

1. X (B y = y (B X. 

2. If (Tj and ag commute, then i^®^ = i^®^; otherwise, 

■S0j? _ _AyQx 



Proof: Without loss of generality, we can assume that 
tr(C/C/-^) is real, because, if otherwise, we can replace U 
by some Ui € [U] such that tr([/if7f ) = \tT{UU'^)\ and 
D{U,On) = D{Ui,On)- 

Let A = {U + U*)/2. Then A is real and 



= ([/t + [/T)(^^^*)/4 

= (2/ + C/tC/* + C/^;7)/4 



(26) 



Moreover, 



\\Ar = tr{AU) 

= ti{2I + WU* + U'^U)/4: (27) 
> N{l-S/2). 

Suppose A has the singular value decomposition A — 
ViJ^V2, where Vi,V2 € Ojv since A is real. Then define 
V — V1V2 G On- Then similar to the proof of lemma [TJ 
we can show 



DiA, V) < VS/2. 



Meanwhile, we have 



\U~A\\^ = ||(L/-C/*)/2||2 

= tr((C/t -C/^)(C/-C/*))/4 

= tT(2I-WU* - [/^t/)/4 

< NS/2, 



which implies 



D{U,A) < VS/2. 



Therefore, 



(28) 



(29) 



(30) 



D{U,On) <D{U,V) 

<D{U,A) + D{A,V) (31) 

< Vs. 



III. TESTING ORTHOGONAL GROUP 

Let us begin with the testing of orthogonal group 

On = {U eVN ■■ U'^U = UU'^ = /} 
= {U eVN -U = U*}, 



(25) 



which is an important subgroup of Ua? that has wide 
applications. Our first result is 

Theorem 2. The orthogonal group On can be e-tested 
with query complexity 0{l/e^). 

Proof: Our basic idea is to show that: (1) if UU'^ is 
close to /, then U is close to Oat; (2) we can test the 
proximity between UU'^ and / very efficiently. Now we 
provide the details. 
First, we prove 

Lemma 3. If U e Ujv satisfies |tr(C/C/^)| > N{1 - 6), 
then D{U,On) < VS- 



Now consider the following testing algorithm: 



Algorithm 1 Testing On 

Input: U £ Ujv is given as a blackbox. e > is a proximity 
parameter. 

Steps: 

1. Repeat the following procedm-e 0(l/e^) times 

• Prepare the state |$^); 

• Apply the operation C/ ® [/ to it; 

• Perform the measurement {|'l'jv)("I'jvli ^ ~ 

on the state. If the outcome cor- 
responds to then this iteration is 
successful. 

2. Accept if and only if all iterations are successful. 
Note that 

(/0A)|$+) - (A^0 /)!$+), VAgMat^w. (32) 
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So 

{U<E>U)\^+) = {UU^^I)\<^+). (33) 
The probability that an iteration is successful is 



Algorithm 2 Testing fc- Juntas 



= |tr(f/[/^)|VA^' 
= (l-7)^ 



(34) 



assuming |tr([/J7^)| = A^(l — 7) for some < 7 < 1. 
Suppose we run iterations in step 1. Then the algo- 
rithm accepts U with probability (1 — 7)^/' < e"'^'^!^ , 
which is smaller than 1/3 for any 7 > e'^ . Hence if the al- 
gorithm accepts accepts V with probability at least 1/3, 
then 7 < and by lemma [3] we have D{U,On) < £• On 
the other hand, for any U £ [Oat], obviously every itera- 
tion is successful and the algorithm always accepts it. So 
this algorithm e-tests On by making 0(l/e^) queries to 
U. m 

Note that — l^'(o'g)) can be efficiently pre- 

pared, and the measurement {|$^)($^|,/- 
can be simulated by a finer measurement in the basis 
{|w(crj))}jgzj. Hence, algorithm [1] can be efficiently im- 
plemented. 



IV. TESTING QUANTUM JUNTAS 

Given a unitary operator U acting on an n-qubit sys- 
tem, we might want to know if it only nontrivially acts on 
at most k of the n qubits. Formally, let [n] = {1, 2, . . . , n} 
be the set of indices of the qubits. For any T C [n], we 
also use T to denote the subsystem composed of qubits 
whose indices are in T. We say that U € Un only non- 
trivially acts on subsystem T ii U — Vr ® It" for some 
V E U2IT1 , where Vr and /t^ indicate that V and / act 
on the subsystem T and T"^ = [n] \ T respectively. The 
set of quantum fc-juntas is defined as 



fc- Juntas = {[/ e Utv : 3r C [n], |T| = fc, s.t. 

U ^Vt <S) It" for some V G Usfc}. 



(35) 



Theorem 4. k-Juntas can he e-tested with query com- 
plexity 0(fc log(fc/e)/e^). 

Proof: Our basic idea is to consider the Pauli decom- 
position given by Eqs. (ll9ll23| . For any a? G Z", let 
supp(x) — {i ■ Xi 0} and \x\ — |supp(a;)|. Note that if 
U — Vt ® It" for some T C [n], then fix{U) = for any 
X with supp(a;) ^ T. So for any fc-junta U, if we mea- 
sure the state \v{U)) in the basis {|u((7j))}2:gzj, then we 
only obtain outcome x satisfying supp(x) C T for some 
T C [n] with \T\ = k. The difficult part is to prove the 
converse is also true in the approximate sense. Namely, 
if we obtain outcome x satisfying the same condition for 
sufficiently high probability, then U is close to a fc-junta. 

Now we give the details. Consider the following testing 
algorithm: 



Input: U G Ujv is given as a blackbox. e > is a 
proximity parameter. 

Steps: 

1. Let W = 0. 

2. Repeat the following procedure 
0(fclog(fe/e)/e^) times: 

• Measurement the state \v{U)) in the ba- 
sis {\v{as))}^^z2- Suppose we get the 
outcome x. 

• Update W W U supp(3;). 

• If |M^| > fc, then reject and quit. 

3. If none of the iterations rejects, then accept. 



Obviously, this algorithm accepts all fc-juntas. It re- 
mains to show that if it accepts U with probability at 
least 1/3, then £)([/, fc- Juntas) < e. To prove this, it is 
enough to prove the following statements: 

Let = {f e Z!^ : supp(f) C T} for any T C [n]. 
Then 

1. if E lA*5(t/)P < l-eV4foranyTC [n], |T| < fc, 

then the algorithm accepts U with probability at 
most 1/3; 

2. if J2 If^siU)]"^ > 1-5 for some T C [n], then there 
exists V e U2IT1 such that D{U, Vr <^ It-) < ^VS. 

The desired result follows immediately from the two 
statements. 

To prove the first statement, consider the following 
classical game. Repeat the following procedure m times: 
each time, we randomly sample a string from ZJ such 
that any a; G Z" is chosen with probability = |/zj([/)p. 
Let xi,X2, ■ ■ ■ ,Xm be the samples. We win the game if 
and only if 



I U supp(fj)| < fc. 
i=i 



If for any T C [71], |T| < fc. 



Pr[supp(.f) C T] = Ps<l-S, 



(36) 



(37) 



then what is the maximal probability of winning the 



game? 

Of course, we can give a simple upper bound (^)(1 — 

m 

(5)™, since there are (^) possibilities of = IJ supp(S'j) 

such that \W\ < fc, and when W is fixed, Pr[supp(ai'j) C 
W] < 1 — (5 for each j. However, this bound is not good 
enough since it depends on n. To get a better bound, 
we introduce the following concept: given a sequence of 
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samples xi,X2, ■ ■ ■ , Xm, we say that j € [m] is a support- 
defining position if 



supp(fj) 2 U supp(fj). 



(38) 



Note that if < k, then there can be at most k 
support-defining positions. Consequently, we can find 

k 

1 < ji < j2 < • • • < jfc < ™ such that = y supp(xjj. 



1=1 



No matter what W is, for any j € J2, ■ • • i Jfc}, we 

always have Pr[supp(a;j) C W] < 1 — 5. Since there can 
be at most (™) choices of ji , j2 , ■ • ■ , jk , by union bound 
we get 



Pv[\W\<k] < (™)(l~(5)™-'= 

^ gfc log m— 5(m— /c) 



(39) 



By choosing m = 0{k\og{k/d)/S), we can make this 
probability smaller than 1/3. Setting S = yields 
statement 1. 

Now we prove the second statement. Let 



^ - E ^^s{U)crs: ^ At ® It-= 
B = U-A= Y: f^AU)<Js. 



(40) 



Then we have 
which implies 



N E \f^siU)\'>Nil-6), 



||i||2 > 2l^l(l-(5). 



(41) 



(42) 
(43) 



Furthermore, let lip) be an arbitrary state on subsystem 
T and let p be the uniformly mixed state on subsystem 
T"^. Then we have 

1 =: tr{WU{\i^){i^\(E)p)) 

= tr(AtA(|V')(V'l P)) + tr(StB(|^)(V'| p)) 

+tr(AtB(|V') (V'l <E)p))+ tr(Bt A(| V') (V-l ® p)) 
= tr(itA(|^)(V'|) + tr(BtS(|V)(V'| ® p)) 
> tr(iti(|^)(^|), 

(44) 

where in the third step we use the fact 

ps{A^ B) ^ ^,s{B^ A) ^ 0, yxeZj. (45) 

Since \ip) is arbitrary, we get 

iti < /. (46) 

By Egs. p^ . P5|) and lemma [1] we get there exists some 
V gV21t\ such that 



and hence 

D{A,Vt It") < \/S/2- 
Meanwhile, we have 

D{U,A) <^\\U-A\\ 

< 7^2- 

As a result, 

D{U, Vt «) ItO < D{U, A) + D{A, Vt ® h-) 

= 2VS. 



(48) 



(49) 



(50) 



Note that algorithm [2] can also be efficiently imple- 
mented. 



V. TESTING CLIFFORD GROUP 

The Pauli group on n qubits is defined as the subgroup 
of Uat generated hyXi=X(g)I(g)---(g)I,Zi^Z<g>I(g) 
•••«)/, X2=/®X(g)---(g)/, Z2=/(g)Zi8)---(8)/, 
X„ = I <E)I(E> ■ ■ -(E) X, Zn ^ 1(g) I(g> ■ ■ ■<S) Z. Equivalently, 



(51) 



The Clifford group on n qubits is defined as the normal- 
izer of Vn, i.e. 



Cn = {U ellN ■■ UhW e Vn, yh e Vn}- 



(52) 



Both Pauli group and Clifford group p lay important roles 
in quantum error correction |32h34| | and fault-tolerant 
quantum computation (ssl . [36j . 

Before stating our result about testing Clifford group, 
it is necessary to first present the following result about 
testing Pauli group. 

Lemma 5 (implicit in Ref . [23| ) ■ The Pauli group Vn can 
be e-tested with query complexity 0(l/e^). 

Proof: Consider the following testing algorithm: 

Algorithm 3 Testing Vn 

Input: U G Ujv is given as a blackbox. e > is a 
proximity parameter. 

Steps: 

Prepare 0(l/e^) copies of \v{U)). Measure 
each copy in the basis {|i'(o"s))}sgzj. If all 
measurements get the same outcome, then ac- 
cept. Otherwise, reject. 

Let m — c/e^ be the number of copies of \v{U)) used 
in the algorithm. Since 



D{A,V) < ^/Sj2, 



(47) 



HU)) = E l^AUMx)), 



(53) 
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the algorithm accepts ?7 with probability ^ 



2m 



E U ^ a So for some xq G ZJ, then l^j,-, = 1 and 
— for any x ^ xq, which implies that every 
measurement gets the outcome xq and the algorithm ac- 
cepts U with certainty. On the other hand, if D{U, Vn) > 
e, then by Eqs.® and ^ we get \^JLx{U)\ < 1 - for 
any x € Z4. Consequently, 

1. ifO<e<l/ \/2 , then the probability of the algo- 
rithm accepting U is 



E lMs(C/)r™ <(l-e2)2™ + (e 

< 2(1 -e2)2m 



^2^2m 



(54) 



By choosing a sufficiently large c, we can make this 
probability smaller than 1/3. 

2. if 1/a/2 < e < 1, then the probability of the algo- 
rithm accepting U is 

E Mc/)P" <(i-e')''"((T^ + i) 

< 2(1 - e2)2m-2 (55) 

which is smaller than 1/3 as long as m > 2. 

Overall, this algorithm e-tcsts Vn with query complexity 
0(l/e2). ■ 

Now we turn to the testing of Clifford group. Note 
that if i7 G Cn, then D{UasU\Vn) = for any x £ ZJ. 
The following lemma shows that the converse is also true 
in the approximate sense. 

Lemma 6. If D{U(J3U\Vn) < 5 for any x G Z", then 
D{U,Cn) < 46. 

Proof: See Appendix A. ■ 

Therefore, in order to test whether U is in [Cn] or far 
from it, it is sufficient to test whether UasW is in [Vn] 
or far from it for every x G Z". Note that ag can be 
generated by ATi, Zi, . . . , X„, Z„. So UagW can be gen- 
erated by UXiW, UZiW,..., UXnU\ UZ^W. For ex- 
ample, a g = XiZ2XnZn G Vn, then we have UgW = 
{UXiW){UZ2W)iUXnU'<){UZnW). This means that 
UcrgW is close to Vn for any x if and only if UXjW 
and UZjW are close to Vn for every j. This suggests 
that we can apply the algorithm [3] to each UXj and 
UZjW , and accepts U if and only if every subtest accepts 
with an appropriate proximity parameter. An algorithm 
based on a similar idea was given in Ref. ,25]. However, 
this approach has the drawback that it needs to execute 
2n subtests on UXjU^^s and UZjU^'s, and consequently 
its query complexity depends on n. Here we present a 
better algorithm whose query complexity only depends 
on the proximity parameter e. 



One basic idea is that we still test the distance between 
U ugU^ and Vn, but we only do this test for a few random 
x's. The key observation is that UVnU^ is a group. This 
group structure ensures that, UagW is close to Vn for 
any x if and only if UagU^ is close to Vn for a sufficiently 
large fraction of x. Specifically, 

Lemma 7. // at least 2/3 fraction of x £ Z4 satis- 
fies D{UasU\Vn) < 6, then for any y G Z4, we have 
D{UagU\Vn) < 2S. 

Proof: We prove this result by using the pigeonhole prin- 
ciple. A a; G Z4 is said to be 6 -good if D{UagW ,Vn) < 5. 
Fix a j/G Z4. Divide Z4 into 4"/ 2 pairs: each pair con- 
sists of X and x Q) y for some x G ZJ . Since at least 2/3 
fraction of Z" is (5-good, at least one pair contains two 
(5-good vectors. Let z and z(Byhe any such a pair. Then 

D{Ua^W,Vn) = DiUasW ■Uiasa^)W,Vn) 
= D{UasW -Uas^ffU^Vn) 

< D{UasW,Vn) + D{Uas^^U\Vn) 

< 25. 

(56) 
■ 

Therefore, we can run the algorithm [3] on only a few ran- 
dom a; G Z", and then estimate the fraction of good x's 
with sufficiently good precision and hence estimate the 
distance between U and C„. Now we give more details. 

Theorem 8. // both U and W can be accessed, then the 
Clifford group Cn can be e-tested with query complexity 

Proof: Consider the following testing algorithm: 

Algorithm 4 Testing C„ 

Input: U,U'' e Uiv are given as blackboxes. e > is 
a proximity parameter. 

Steps: 

1. Repeat the following procedure 0(1) times: 

• Pick a. X £ Z4 uniformly at random. 

• Run the algorithm|3]on UagU^ with the 
proximity parameter e/8, where the or- 
acle U agU^ is simulated by concatenat- 
ing U\ a 3 and U. 

2. If all iterations accept, then accept. Other- 
wise, reject. 

If C/ G [Cn\, then UagW G [Vn] for any x and hence 
the algorithm always accepts U. On the other hand, we 
will show that if the algorithm accepts U with probabil- 
ity greater than 1/3, then at least 2/3 fraction of x G Z" 
is e/8-good, and then by lemma IH] and lemma [71 we get 
D{U,Cn) < £• Suppose the fraction of e/8-good x is p. 
Then if we pick such a good x, then algorithm 3 accepts 
UagU^ with probability at most 1; otherwise, the algo- 
rithm 3 accepts UagW with probability at most 1/3. 
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Hence, the probability that algorithm |4] accepts U is at 
most (p + (1 — p)/i)'~^ , where C = 0{1) is the number 
of iterations. If p < 2/3, then this probability is at most 
(7/9)'^ which is smaller than 1/3 as long as C > 5. ■ 

Note that algorithm [3] is efficiently implementable, and 
hence so is algorithm 21 



Without loss of generality, we can assume that {U,Wi) 
is real for all i, because, if otherwise, we can replace Wi 
by some W- £ [Wi] so that this condition is fulfilled. 
Then we have 



= 1-D^{U,W^) 



(59) 



VI. TESTING ANY FINITE SUBSET 

So far we have studied the testing of several special 
subsets of Utv- In this section, we will present an algo- 
rithm that tests any finite subset of U^v and thus also 
give an upper bound on the query complexity. 



Theorem 9. Suppose S = {VKi, VFz, . . . , Wm} cVn- If 
N(l — S), then S can be e-tested 

logM 



rnax |(m,W^, 

l<i<j<M 



with query complexity 0{ 



minje^, 5} 



Proof: Our basic idea is as follows. Suppose U is ei- 
ther in [S] or far from [S]. Then \v{U)) is either in 
the subspace spanned by |u(W2)), . . . , \v{Wm)), 

or has bounded projection onto this subspace. Since 
|(<^®/c|^®fc^| = K^'IV')!'', we would expect that 
\v{U))®'' is either in the subspace spanned by 
\v{Wx))®'' ,\v{W2))'^'' , ■ . ■ ,\v{Wm))'^'' , or has exponen- 
tially small projection onto this subspace. So by making 
a projective measurement on \v{U))^^ we should be able 
to distinguish the two cases. In the next, we are going to 
make this argument rigorous. 

Consider the following testing algorithm: 



Moreover, without loss of generality, we can also as- 
sume that [W^i], [W2], ■ • ■ , [W'm] are disjoint. Then 
\v{Wi)),\v{W2))-,- ■ ■ ■,\v{Wm)) are linearly indepen- 
dent and Wk is Af-dimensional. Let = 

M 

^ \,j\'v{Wj))®^}f^^ be an arbitrary orthonormal ba- 
sis for Wk- Then we have 



M 

M = Y.{^^\^^) 
i=l 
M 

= E K,jKr{{^iw,)\v{Wr))Y 

M M M 

= E E |A.,f + E E xi,x^,r{{v{w,)\v{w,,)})^ 
1=1 j=i 1=1 j^j' 

M M M |, |2,|, ,2 



I.e. 



> EEIA.,,P-E E 

i=ij=i 1=1 j^ji 

M AI 

> iEEIA.,f, 

i=ij=i 



M AI 

E E 

i=ij=i 



E E |A.,,f <2M. 



(60) 



(61) 



Algorithm 5 Testing S = {Wi,Wi, Wm} C Vn 

Input: U £ Ujv is given as a blackbox. e > is a 
proximity parameter. And suppose 5 = 1 — 



Steps: 



i max \{Wi,Wj) 

" l<i<j<M 



j^rrjy) copies of the 



1. Prepare K = 
state \v{U)). 

2. Perform the measurement {11/^,1 — TIk} 
on \v{U))^'^ , where Ilif is the projec- 
tion operator onto the subspace Wk = 
span{|i;(W^i))®^^ |w(VF2)>®^, ■ • • , 

If the outcome corresponds to Ilif, then ac- 
cept; otherwise, reject. 

Obviously, HUE [S], then the algorithm always ac- 
cepts. So it remains to show if D{U,S) > e, then the 
algorithm accepts with probability at most 1/3. 

Choose K = 0( '"ft^, ) such that 



K 



(1 



< e 



< e 



-SK 



< 



< 



5M' 
1 

5M' 



(57) 
(58) 



Now if we perform the measurement {IIk,I — Hx} on 
\v{U))'^^ , then the probability of obtaining the outcome 
corresponding to Hk is 



tTiUK{\viU)){viU)\m 

M 
i=l 

M M 

= i:\T.KA{v{w,)\v{u))) 
1=1 j=i 

M M 

= E E xi,\r{Hw,)\v{uW{{v{u)\v{w,))) 

M M 

< (1-O'^E E \Kj\\K.r\ 

M Af ,, |2.|i ,2 
- 2517^ 2 



\K\2 



K 



AI M 



< 2^ E E \XrJ? 



1=1 3 = 1 



< 



10' 



(62) 

where in the fourth step we use Ea.(|59p. and in the last 
step we use Eg. lpTj) . ■ 
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A. Example: Testing Permutations 

Let us demonstrate an application of theorem [9] to the 
permutation group. These operators just relabel the n 
qubits of the system but do nothing else. Formally, let 
Sn denote the group of permutations over {1, 2, . . . , n}. 
Then any r G is viewed as a unitary operator as 
follows 

T\ii,i2,...,in) = |V(1),V(2), ■ • ■ ,V(n)), (63) 

for any ii, 12, . . . , i„ = 0, 1. 

For any n ^ T2, let j — t^^T2 ^ I. Note that 7 can 
be decomposed into several disjoint cycles 

7 — (ai.ii • ■ • I ai,ki){a2,i, • ■ • , 0-2^) ■ ■ ■ • ■ • j o.i,ki), 

(64) 

where I < n — 1. So 

(ti,T2) = tr(T/T2) 

= tr(7) 

= tr( J2 N7(l):---:«7(n))(«l, •■•,««!) 

n,...,i„=04 

= J2 ■ • ■ :*n|i7(l), ■ • ■ ,«7(n))- 

ii,...,i„=0,l 

(65) 

By Eq. ljM]) . the only nonzero terms on the right-hand 
side are those satisfying 

«a,,i = ia,,2 = • • • = ia,_,^ , Vj = 1, 2, . . . , L (66) 

Hence 

(ri,r2)=2'<|. (67) 

This holds for any ri 7^ r2. Therefore, 

<5=l-lmax|(Ti,r2)| > ^. (68) 

Besides, note that |S'„| = n! = 0(e"'°s"). So by theorem 
m we get 

Theorem 10. The permutation group Sn can be e-tested 
with query complexity 0{nlogn/e^). 

Note that it is unknown whether algorithm [S] can be 
implemented efficiently in general. So it still remains 
open to find a testing algorithm for 5„ that has both 
polynomial query complexity and polynomial computa- 
tional complexity. 

VII. CONCLUSION 

In summary, we have systematically studied property 
testing of unitary operators with respect to the distance 
measure D that reflects the average difference between 
unitary operators. We present efhcient algorithms for 
testing the orthogonal group, quantum juntas and Clif- 
ford group. All these algorithms have only one-sided er- 
ror and their query complexities are independent of the 



system's size. We also give an algorithm that tests any 
finite subset of the unitary group, and show an applica- 
tion of this algorithm to the permutation group. This 
algorithm also has one-sided error and polynomial query 
complexity, although it is unknown whether it can be 
efhciently implemented in general. 

Despite the progress made in recent papers and ours, 
the testing of quantum objects (states or operations) still 
remains widely open. We hope that our work can shed 
some light on this topic and stimulate further research. 
Here are several directions that seems particularly inter- 
esting to us: 

First, in this paper we focus on giving upper bounds 
on the query complexity of testing unitary operators. It 
is worth developing powerful techniques that can derive 
lower bounds on the query complexity of the same task. 
In particular, it is interesting to prove that our testing 
algorithms are optimal, or give better algorithms. 

Second, as mentioned above, we do not know whether 
algorithm [5] can be efhciently implemented, or whether a 
better algorithm can be given for testing general discrete 
subsets. And it would be interesting to give an efficient 
algorithm for testing the permutation group. Further- 
more, can we give some general results on testing contin- 
uous subgroups? 

Third, here we only considered testing unitary oper- 
ators. It is worth exploring property testing of general 
quantum states and quantum operations. For example, 
is it possible to test (or estimate) the entangling power 
of multipartite quantum operations by using only a few 
queries? Is it possible to test whether a given multipartite 
state belongs to an interesting class, such as symmetric 
or antisymmetric states, by consuming only a few copies? 

At last, the query complexity of property testing cru- 
cially depends on the distance measure used. Here we 
considered a particular distance measure D that seems 
quite suitable for comparing unitary operators. It would 
be interesting to study the testing of quantum states or 
channels with respect to other natural distance measures. 
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Appendix A: Proof of lemma [5] 

If (5 > 1/4, then obviously £>([/, C„) < 1 < 4(5. So from 
now on we assume 5 < 1/4. Define : Z4 — > Z4 and 
e : ZJ^ -> [0, 27r) such that for any f g 

DiUagUWn) = -^WUasW - e^^'^^^ap^s^W < S, 

(Al) 

i.e. e*®(^Vi;'(x) is the closest element to UagU^ in \Pn]- 
We will prove that 

1. Q{x) = or TT, Vf e I'l- 
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2. F{x)®F{y) = F{x®y),yx,y€Z^; 



Now 



3. as and ag commute (or anticommute) if and only 
if (Tf{x) and (JF(y) commute (or anticommute), 
V£,y G ZJ. 

Let us first assume these statements are true. Then the 
mapping —> 'ysaF{x) is an isomorphism, where 7^ — 
gieix) _ _|_-[^ Hence there exists a Clifford operator C € 
Cn such that 



(All) 



Then by Ea. (|XT|) . we have 

{UasU^CagC'f) > N{1 - S^), Vf e ZJ. 
Consequently, 

< i7^tr((7( 2 asWCas)C^) 

= ^tr(t/(mr(C/tC)/)Ct) 
= j,|tr({/tC)p, 

where in the fourth step we use the fact 

^ J2 TxAas = tT{A)I, VA e 



(A2) 
(A3) 



(A4) 



Hence, we have 



|(C/,C)| >NVI^>N{1-5'^). 
Then by Eq.® we obtain 

D{U, C) < 5. 



(A5) 



(A6) 



(A7) 



Now we prove statements 1-3. The first one is easy. 
Since as is Hermitian, UagU^ is also Hermitian. As a re- 
sult, tr(?7crjC/V^(^)) is real. 9(5;) is chosen from [0, 27r) 
such that 

e-B(-)iy(UasU^aF(s)) = \iv{UasW ap(^3))\. (A8) 

So e*®*^^-* can only be ±1 and Q{x) can only be or tt. 
To prove the second statement, let us suppose 



z^F(x) 
z^Fiv) 



(A9) 



By D{UasU\aF{x)) ^ ^ and D{UajjU\aF(ij)) < S, we 
get 



\^Fiy)\>l-S'- 



for some rj/s. Since D{Uas(^jjU\aF(SBy)) — '^^ have 
hF(Jey)l > (A12) 
Meanwhile, the coefhcient T]F(x)(SF(y) satisfies 

\VF{S)(BFiy)\ = I E «^'''®^'''Aj=^(i-i)'^F(i-2)l 

zi(Sz2=F{x)eF{y) 

+ E i'''^'^^^'^^^^^^^^>^z^mFixnFiyl\ 
z^F{x) 

> I^F(x)CF(a)l - E \^z£,z(BF(x)eF{y)\ 

z^F{x) 

> _ jj2^2 _ ^ |Ai-|^+|C?e-F(s)e-F(ff)l^ 



> (1 - J2)2 _ (1 _ ^2)2) 2 

> 1-4(5^. 



If F{x) ® i^(y) 7^ F{x ® j7), then we would have 
1 > l?7F(xej?)P + l?7F(x)eF(a)P 

> (1 - (52)2 + (1 „ 4^2)2 

> 2 - 10,52 

> 1, 



(A13) 



(A14) 



which is a contradiction. Therefore, F{x)®F{y) = F{x(B 

y)- , 

Finally, we prove statement 3. Let us first consider the 
case when a^ and ag commute. Then UagU^ and Uagll^ 
also commute. We still assume Ea. (|A9p . and suppose 

UasW-Uay-W = E W 



(A15) 



(AlO) 



If crF(K) and cri?(j;) anticommute, then we have 
•F(£)oF(y) ^ _jF(i?)0F(J)^ and furthermore, by Eq.(|Xl3|, 
we get 

= \VFix)<SF{y) - XF(y)(!)F(x)\ 

^ \^,FimF(y).,Fly)eFlx)^XFix)^^iy) 

+ E (z^®(^®-'^(^)®-'^(^)) - i(^"®-P(^)®-P'(i?))0-?) 
•^i'Czef'(5)eF(i?) I 

> '^\^F{x)S,F(y)\ ~ E l''^z'?zeF(a)ef'(i/)l 

z^F{x) 

> 2(l-4<52) 

> 0, 

(A16) 

which is a contradiction. So aF{x) and aF{y} must com- 
mute. The case when as and aij anticommute can be 
handled by a similar argument. 
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